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Thermoelectric Cooler Design is a Microsoft Windows program to be used as an aid in 
the design of thermoelectric cooler devices. This program was written to be used to 
quickly model and compare alternative designs. A couple's optimum coefficient of 
performance and maximum heat pumping can be quickly determined. Other major 
features of the program include the ability to change material properties and dimensions 
of couples, analyze cascaded couples, and graph performance parameters. A brief 
description of thermoelectric cooler theory, modelling assumptions and complete source 
code listing is included. 
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I. INTRODUCTION 



A. FEATURES OF THERMOELECTRIC COOLERS 



The unique features of thermoelectric coolers include: 



1 . Compact. 

2. No moving parts. 

3. The coefficient of performance is independent of capacity. 

4. The cooling rate is easily controlled. 

5. Reversible operation is possible allowing heating as well as cooling. 

6. Vibration-free operation. 

7. Precise temperature control. 

B. SPACEBORNE APPLICATIONS 

Thermoelectric coolers offer an alternative solution to cooling and heating 
requirements of spacecraft. Spacebome applications for thermoelectric coolers include: 



1. Equipment cooling -- removal of waste heat in order to maintain 
components at their optimal operating temperature. Thermoelectric coolers 
are especially useful for the control of small hot spots. 
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2. Life support temperature control — thermoelectric coolers are very reliable 
making them ideally suited for long duration missions such as the space 
station or manned Mars mission. 

3. Cooling of infrared detectors -- thermoelectric coolers eliminate the need 
for coolant and associated support equipment. Of even more importance is 
the elimination of consumable coolants which also increases the life of the 
spacecraft. 

4. Microgravity applications -- thermoelectric coolers offer precise, 
vibration-free temperature control for free-flying orbiting microgravity 
research platforms. 

Reliability, the lack of moving parts or consumables, and vibration-free operation 
are especially important to spacebome cooling systems. 
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II. THERMOELECTRIC COOLER THEORY 



A. THERMOELECTRIC COUPLE CONFIGURATION 

In it's most basic form, a thermoelectric couple consists of two dissimilar elements 
with contact at two junctions. Figure 1 shows a simple thermoelectric circuit consisting 
of two elements, physically joined on one side while electrically joined at the other. 




Side 



Figure 1 Thermoelectric Cooler 
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B. THERMOELECTRIC EFFECTS 



The background for this section was gathered from Kraus [Ref. 12:pp. 294-303] and 
Ioffe [Ref. 2:pp. 96-128]. Three thermoelectric effects are known; the Seebeck, Peltier, 
and Thomson effects. 

1. Seebeck Effect 

Consider the junction of two materials, A and B. If a temperature differential 
exists between the two junctions, one junction may be referred to as the hot side and the 
other the cold side. An open circuit voltage can be measured between the two junctions. 
This net conversion of thermal energy into electrical energy under zero current 
conditions is known as the Seebeck effect. The voltage observed, referred to as the 
Seebeck voltage, is proportional to the temperature differential between the hot and cold 
sides of the couple: 



dEs « dT 



dE s = ± a dT 




( 1 ) 



where 



Es = Seebeck voltage (V) 
a = Seebeck coefficient of material (V/°C) 

In the case of the couple composed of materials A and B, 
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where 



a,j = Seebeck coefficient of material A 
(Xg = Seebeck coefficient of material B 

2. Peltier Effect 

The Peltier effect is localized at the junctions between materials A and B. When 
a current flows through the junction, heat is either generated or absorbed depending on 
the direction of the current flow. The amount of heat is proportional to the current and is 
known as the Peltier heat: 

dQp Idt 



dQp = ± nldt 



qp = ±nl (W) (2) 

where 

Qp = Peltier heat (J) 
qp - Peltier heat flow (W) 

I = Current (A) 

n = Peltier coefficient of material (V) 

3. Thomson Effect 

This effect concerns the absorption or generation of heat in a conductor carrying 
a current in the presence of a temperature gradient. The amount of heat is proportional to 
the current and temperature gradient: 

dQ t °c I dTdt 
dQ t = ±al dTdt 
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dq t = taldT 



q t = ±C/ J dT (W) 



0 ) 



where 



Q t = Thomson heat (J) 
q t = Thomson heat (W) 

/ = Current (A) 

a = Thomson coefficient of material (V/ C) 



Because the product of o and AT is a voltage, the Thomson voltage is therefore: 



T 2 

E( = ± J cc/r (V) (4) 

T \ 

C. IRREVERSIBLE EFFECTS 

1. Joule Effect 

The evolution of heat when a current passes through a material may be referred 
to as the Joule heat flow: 

qj = I 2 R (W) (5) 



6 



In the case of a couple, the elements A and B are connected in electrical series 



regardless of the direction of the current flow: 

( 6 ) 

(£2 cm) 

(cm 2 ) 

2. Fourier Effect 

The conduction of heat through a material is given by: 



R = PA±l + PB±B_ (n) 
a a a b 

where 

Pj , pg = electrical resistivity of materials A and B 
La, Lb = Length of materials A and B (cm) 

A a, a b ~ Cross sectional area of materials^ and B 



q f = T (W) 



(7) 



and, in the case of a couple with the elements in thermal parallel, 



k A A A k B A B 
^B 



(W/°C) 



where 

k A' k B = Thermal conductivity of materials A and B (W/cm °C) 
La, Lb - Length of materials A and B (cm) 

A A' A B - Cross sectional area of materials /I and B (cm^ ) 



( 8 ) 
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D. NET HEAT ABSORBED AT COLD JUNCTION OF COUPLE 



In this analysis, the Thomson heat will be neglected which greatly simplifies couple 

design procedure without any significant error Kraus [Ref. l:p. 300] and Cadoff [Ref. 

3:p. 8], To determine the net heat pumped by the couple, the temperature distribution 

along the elements of the couple must be determined. Consider a single thermoelectric 
element as indicated in Figure 2, where Ax is the cross sectional area of the element 
normal to the flow of current and heat. 
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Here, the boundary conditions are: 



T = To at x = 0 
T = Te at x = L 



Assuming steady state temperature conditions, there are three heat quantities which 
pertain to an energy balance over the incremental length. Ax: 

1 . The Fourier heat entering at point x: 




2. The Joulean heat generated in the increment Ax : 



3. The Fourier heat leaving at point x + Ax : 

1x+Ax = - kAx [&] x+Ax 



The energy balance in the steady state requires that the heat leaving must be equal to 
the sum of the heat entering and the heat generated in the increment Ax : 



Qx + (}j - q x+Ax 

Substituting for qx, q x +Ax anc * qj gives 
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-“■[«],• 4ir - -"■[«] 

or with some algebra, 

kAx 



(\dL\ 


i 


\l dx Jx+Ax 


IdxJx) 



x+Ax 



I - p Ax 

+ r- = 0 



Ax 



Dividing through by Ax yields 



kA x 



T — 1 


W 


L dx Jx+Ax 


IdxJx 



Ax 



1 1 p 

+ ~~A = 0 

Ax 



And taking the limit as Ax — > 0 gives 



limit 
x — > 0 



1 1 

X 


y-8 
* « 
V. 


ix+Ax 




X 

1 1 


1 

Q_ 




Ax 


) 


Ax 



= 0 



which can be adjusted to 



kA x + 4^ = ° 
dx 2 A x 



(9) 



The general solution to this simple differential equation is obtained through the use of 
a double integration: 

d 2 T , I 1 P* 



dx 2 



kA 



2 + C 1 - 0 
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and 



I 2 p * 2 

T + + Cj* + C 2 = 0 

2k AZ 



Applying the boundary conditions for the hot side where x = 0 and T = To'. 

C 2 = - T 0 



and 



/ 2 p x 2 

T- To + P T +0^ = 0 

2kA x 



At the cold side of the element where x = L and T - Te '. 



I 2 pL 2 

Te - To h ~ — i- C jL = 0 

UAi 



The solution for C j , the first constant of integration, is seen to be 



r _ To- T e I z pL 
1= L 



so that the particular solution of this simple, second order, differential equation which 
gives the temperature profile as a function of distance along the element: 



/ 2 p AC 2 

T - To + — V — T + 

UA 2 



(To~Te\ x 


( o \ 

I 2 pL 


l L r 


K 2.kA x ) 



x = 0 



or 
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T(x) = - 



,2 \ 
/ z p 


J2 


' To -T e 


/ 2 p L ^ 


^2 kA 2 j 


x - 


L 

\ 


2kAl) 



x + To 



( 10 ) 



which is seen to be parabolic. 



Figure 3 illustrates a temperature distribution for a typical element with the following 
parameters: 



z = 3.08 x 10 -3 (°C 1 ) 


(id 


R = 1.73 x 10 -3 (Q) 


(12) 


K = 3.64 x 10 -2 (W/°C) 


(13) 


a = 4.4x10 ^ (V/°C) 


(14) 


T c = 308.6 K 


(15) 


> 

H 

II 

o 

0 

n 


(16) 



Due to the parabolic temperature distribution, the Joulean heat will be transferred to both 
the cold and hot junctions of the couple. 

Differentiating equation (10), setting the derivative equal to zero, 



dT 

dx 



I 2 pX _ To -T e + / 2 P L 
kA 2 1 2kAl 



permits a solution for x = X which is the position where the maximum temperature 



occurs: 



X = 



/ 2 p L 

2kA 2 x 



To - Te 

L 



fipX_ 
kA. I 



12 



360 




X a-iniejaduiaj, 
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Rearranging and substituting A T for To - Te : 



L k_ATA$ 

2 7 2 pZ, 

Assuming x = L and defining / as the fraction of the Joulean heat transferred to the 
cold junction: 

f _ L-X 1 kJXT_A ij 

J L 2 / 2 pZ 2 

or 

<i7) 

where K and R are defined by equations (6) and (8), respectively. 

The net heat absorbed at the cold junction of the couple is equal to the summation of 
the Peltier heat, minus the fraction of Joulean heat transferred to the cold junction: 

9 net = a’Tc - ffi R 



and upon substitution of the expression for / given by equation (17): 



q net = a ITc ~ \i 2 R ~ KAT (W) (18) 



Therefore the net heat pumped is the Peltier heat developed at the cold junction less 
the sum of one half the Joulean heat produced in the couples and the heat due to 
conduction from the hot junction to the cold junction. 
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III. THERMOELECTRIC COOLER MODEL 



A. INTRODUCTION 

The equations in this chapter were used to model the thermoelectric cooler couple 
operating parameters and were given by Kraus [Ref. l:pp. 303 - 326]. 



B. HEAT PUMPING 



1. Current Yielding Maximum Heat Pumped 

The equation for net heat pumped 



q net = a ITc -jI 2 R- KAT (W) 



( 18 ) 



may be differentiated with respect to current, I. When the derivative is set equal to zero, 
the current that yields maximum heat pumped is obtained: 



= oc 7c - IR = 0 

dl 
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Solving for /: 



a T c , A . 
hn = —r— (A) 



( 19 ) 



2. Maximum Heat Pumped 



Substituting Ini into the equation for the net heat pumped yields the maximum 
heat pumped: 



qm 



« 2 

2R 



- KAT (W) 



( 20 ) 



C. MAXIMUM TEMPERATURE DIFFERENTIAL 

Heat pumped as a function of current is shown in Figure 4 using the parameters 

defined in equations (11) through (16). Setting equation (20) which is the maximum heat 
pumped equal to zero and solving for AT, one obtains the maximum temperature 
differential between the hot and cold sides of the couple: 



AT m 



a 2 Tl 

2KR 



(K) 



( 21 ) 



D. OPTIMUM RATIO OF ELEMENT CROSS SECTIONAL AREAS 

The maximum temperature differential given by equation (21), can be increased 
by minimizing the product KR from equations (6) and (8): 
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o 

CD 




a[dnoo/( M ) b 
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Figure 4 Heat Pumped as a Function of Current 



KR 



P A 1 A + P B l _B ' 
A A A B ) 



r j ^A 

L A 



L B 



and if L A = L B an expansion gives 



K R = k A p A + k B p A 




+ k A p B 




+ k B P B 



( 22 ) 



(23) 



An optimization of equation (23) can be accomplished by finding the point 
where the derivative with respect to the area ratio, A A /A B , vanishes: 



dKR 

d{A A !A b ) 



= k APB ~ k BPA 



K A B j 



= 0 



(24) 



Then, a solution of equation (24) yields the area ratio that minimizes KR. Thus, in turn 
yields the optimal ATm '■ 



£i = PA k B 
A B V p B k A 



(25) 



The area ratio given by equation (25) when substituted into equation (22), yields 
the optimum value of KR, designated by <|> 



or 



$ - k A Pa + 2 iWS + k B p B 





(26) 
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Equation (16) has the restrictions that = Lg and that the area ratio be determined 
using equation (25). 



E. FIGURE OF MERIT OF ELEMENT MATERIALS 

Substituting <)> , equation (26), into the equation for A Tm equation (21) : 



or 



ATm = ? \ t2 c 



A T m = 



2 Zlc 



where z, the figure of merit of the materials, is defined as: 



z = 



a" 



<t> 



which may be rewritten as: 



z = 



or 



( j k APA + J k B P B ) 



When = Lg and the area ratio determined from equation (25), the figure of 
merit is a function of the properties of the elements. To maximize the figure of merit, the 

materials used should have a high Seebeck coefficient, low thermal conductivity, and low 

electrical resistivity Cadoff [Ref. 3: p. 21]. 
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